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Abstract 

Starting from the tri-Hamiltonian formulation of the Lagrange top in a six-dimensional 
phase space, we discuss the possible reductions of the Poisson tensors, the vector 
field and its Hamiltonian functions on a four-dimensional space. We show that the 
vector field of the Lagrange top possesses, on the reduced phase space, a quasi- 
bi-Hamiltonian formulation, which provides a set of separation variables for the 
corresponding Hamilton-Jacobi equation. 
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1 Introduction 



The classical theory of separation of variables for the Hamilton- Jacobi equation provides 
the most effective tool to solve the equations of motion of a given Hamiltonian system. 
In this framework, the main problem is to have an efficient (possibly algorithmic) way 
to produce a set of separation variables. To this purpose, two new approaches, stemming 
from soliton theory, have been recently introduced: the "magic Sklyanin's recipe" [|1|, based 
on the Lax representation of the equations of the motion, and the bi-Hamiltonian (bH) 
approach to separation of variables ^ ^, |^, based on the bi-Hamiltonian structures 
associated with the equations of motion. A remarkable feature of the latter approach is 
that if the Hamiltonian system admits a quasi- bi-Hamiltonian (qbH) formulation, then 
a set of separation variables can be algorithmically computed 0; moreover, the qbH 
property is independent of the coordinate system in which the bH structure is written 
down. 

The aim of this paper is to apply the approach based on the qbH property to the classical 
Lagrange top (LT); in particular, we show how the (complex) separation variables for LT, 
introduced in in an algebraic-geometric setting, arise quite naturally as distinguished 
functions for its tri-Hamiltonian structure. 

The starting point of our analysis is the fact that, on a six-dimensional phase space M, 
the LT vector field admits a tri-Hamiltonian formulation = P^dha (throughout 
the paper, the index a takes values 0,1,2), each one of the three compatible Poisson 
tensors Pa possessing two independent Casimir functions. 

When one tries to eliminate the Casimirs by fixing their values, one is faced with a typical 
situation, occurring also for other bH finite-dimensional integrable systems |0, |], to 
each one of the symplectic leaves Sa one can restrict only the vector field Xl and the 
corresponding pair (P^, h^), but not the entire triple of the Poisson tensors, so that the 
tri-Hamiltonian formulation of Xl is lost under restriction. Nevertheless, using a more 
general reduction process a la Marsden-Ratiu, we will show that the symplectic leaf 5*0 of 
the Poisson tensor Pq can be endowed with a Poisson- Nijhenuis structure P, [l^ (hence a 
bH structure) and that X^ can be given a qbH formulation. So, the separability of LT is 
obtained from its Hamiltonian structures as a natural outcome of the reduction process. 
The paper is organised as follows. In Section ^ the tri-Hamiltonian structure of LT is 
shortly reviewed; in Section ^ the main properties of the qbH model are discussed in 
view of application to LT. In Sections ^ and ^, respectively, the reduction of the Poisson 
tensors P^ and of the vector field Xl with its Hamiltonian functions are considered; 
the qbH formulation for Xl is explicitly constructed, together with a solution of the 
corresponding Hamilton- Jacobi equation. Our results are summarised in Section 6, where 
some potential extensions of this work are pointed out. 
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The multi-Hamiltonian structure of the Lagrange 
top 



A modem formulation of LT can be found in [ITTl |T2i: as usual in this framework, the com- 



ponents of vectors and covectors and the entries of matrices are referred to the comoving 
frame, whose axes are the principal inertia axes of the top, with fixed point O. 
The phase space M of LT is parametrised by the pair m = (w, 7), where u = {ui,U2,uj3)'^ 
and 7 = (71,72,73)^ are the angular velocity and the vertical unit vector, respectively. 
The following notations are introduced: fi is the mass of the top, g the acceleration of 
gravity, J = dia.g{A, A, cA) the principal inertia matrix (c 7^ 1), G = (0, 0, a)^ is the 
center of mass; at last, normalisations are chosen so that fiag/A = 1. 
The Euler-Poisson equations are dLo/dt = Mo (change of the angular momentum) and 
d'~f/dt = (invariance of the vertical unit vector); with the above notations and normali- 
sations, these equations take the well-known form 



dm 
~dt 



Xdm) 



Xdm) 



( (1 - C)u)2^'i - 72 \ 

-(1 - c)u;3u;i + 7i 




V 



72^-^3 - 73^2 

73u;i - 7iu;3 

7lCo'2 — 72t^l 



(2.1) 



The LT vector field Xi^ can be given a tri-Hamiltonian formulation 

= Podho = Pidhi = P2dh2 ; 
the compatible Poisson tensors Pa, written in matrix block form, are 



B 
B C 



Pi 



-B 
F 



P2 



T R 
-R^ 



(2.2) 



(2.3) 



where B, C, T, T and R are 3x3 matrices 

c UJ3 
-c uj^ 

UJ2 ~UJl 
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—UJ2 
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C LU^ 
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-UJ2/C UJi/c 

The Hamiltonian functions can be written as 
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(2.4) 



ho = -F4 + 2acFiF3, 



F3 - 2acFiF2 



(2.5) 
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where o = ^ and 

Fi = ^3, F2 = ^{ul + ul + cul)-js, (2.6) 

F3 = Uiji + UJ2J2 + C , -^4 = 7l + 72 + 73 • 

As it is known, the functions Fi {i = 1,...,4) are integrals of motion for Eq.( |2.1| ); they 
are independent and in involution w.r.t. each one of the three Poisson tensors. Moreover, 
(^1,^2) are Casimir functions of Pq, (^1,^4) of Pi and (^3,^4) of P2. 
The vector field can be immersed in two different bH chains, starting and ending with 
the Casimirs of the Poisson tensors Pa- 

PodF2 = 0, P2dF2 = Podho = Xl, P2dho = Pod{-aF^), P2d{-aF.^) = ; (2.7) 

PodF2 = 0, PidF2 = Podhi, P^dhi = Podho = Xl, Pidho = Pod{-acFiFi), 

Pid{-acFiFi) = . 



2.1 Remark. The Hamiltonian formulation of LT w.r.t. P2 is classical (see, e.g.,^3^]). 
The hH formulation w.r.t. (P07-P2) was introduced in / [7^ / in the semidirect product 
so(3) X so(3), and was later recovered in ^ in an algebraic- geometric setting. The tri- 
Hamiltonian formulation w.r.t. (Pq, -Pi, P2) was constructed in by a suitable reduction 
of the Lie-Poisson pencil defined in the direct sum of three copies o/so(3). (To compare 
the above-quoted results, let us recall that the angular momentum and the vertical unit 
vector are taken as dynamical variables in / (7^ , \TB^ , [7^ , whereas the angular momentum 
is replaced by the angular velocity uj in and in the present paper.) o 



3 The quasi-bi-Hamiltonian model 

The qbH model was introduced in [T^, || and developed in |P, [T^ (see also and references 



therein). Here we summarise some facts to be used in the rest of the paper. 

Let Qo, Qi be two compatible Poisson tensors on a manifold M; a vector field X is said 

to admit a qbH formulation w.r.t. Qq and Qi if there are three functions p, H, K such 

that 

X = Qo dH = - Qi dK . (3.1) 
P 

In other words, X is Hamiltonian w.r.t. Qq with Hamiltonian function H, and it is quasi- 
Hamiltonian (qH) w.r.t. Qi, with qH function K and conformal factor 1/p. In spite 
of the presence of p, equation ( p.l|) implies that H and K are in involution w.r.t. both 
Poisson brackets corresponding to Qo and Qi (as well as in the particular bH case p = 1). 
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If dim M = 2n, the qbH formulation is said to be of maximal rank if at each point m G M 
the Poisson tensors Qo, Qi are non degenerate and the associated tensor N = Qi Q^^ 
(with vanishing Nijenhuis torsion) has n independent eigenvalues Ai(m), A„(m). In this 
case, one can introduce a local chart {Xi, ^i) {i = 1,2, ...,n), called a Darboux-Nijenhuis 
chart [0, such that Qq, Qi and take the canonical form 

--{ti)- (-) 

with A = diag(Ai, A^); in general, the coordinate functions /ij, canonically conjugated 
to Ai, can be computed by quadratures. At last, the qbH formulation is said to be of 
Pfaffian type if p = HILi -^i- 

The following result has been proved in for a Pfaffian qbH vector field. 

3.1 Proposition. The general solution of Eg. ]\ ) for the Pfaffian case is given by 
functions H and K which, in a Darboux-Nijenhuis chart {Xi,fii), take the "canonical" 
form 

H = J2h ^ = Evi' A. = n(A.-A,), (3.3) 

1=1 * i=i « « j^i 

where each fi is an arbitrary function, depending at most on the pair {Xi,fii). Moreover, 
the Hamilton- J acobi equations for both H and K are separable. o 



This Proposition has a straightforward consequence. 

3.2 Corollary. Let X = Qq dH be a Hamiltonian vector field; if in a Q^-Darboux chart 
(x, y) the Hamiltonian H takes the canonical form l \3.!^ , then X admits a Pfaffian qbH 
formulation w.r.t. a Poisson tensor Qi and a qH function K of the form (\3. ^ ) and ( \3. Sj ), 
respectively. 

Viceversa, let X = 1/p Qi dK be a qH vector field w.r.t. Qi; if, in a chart {x,y), Qi and 
K take the canonical forms ( \3.2[ ) ^3.!^ and p = Y[i=i then it is also X = Qq dH with 
Qq and H given by ^3.2J (^. Sj ), respectively. Hence, the chart (x, y) is a Darboux-Nijenhuis 
chart for the Poisson pair Qq, Qi. o 

For n = 2, this Corollary can be slightly generalised, in a way that is useful for subsequent 
applications to LT . 

3.3 Proposition. Let S be a four- dimensional manifold and Y = Qq dH be a Hamil- 
tonian vector field w.r.t. a non degenerate Poisson tensor Qq. Let there is a Darboux 
chart (x, y) such that the Hamiltonian H can be written as a linear combination of two 
functions H, K with the canonical form ( ^. 3^ ), i.e., 

H{x,y) = (3H{x,y) + k{x,y) (3 = const , (3.4) 
H{x,y) = ^ (fi{xi,yi) -/2(x2,y2) 
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K{x,y) = — ^ — (x2fi{xi,yi) -Xif2{x2,y2) 

Xi - X2 V 

Then, the vector field Y admits the Pfaffian qhH formulation ^3. 1\ )- ^73[ ); a Darhoux- 
Nijenhuis chart (A, /i) is given by the following map: 

$ : h-. (A,/i) \ = ^—, = -y^(x, + (^ = 1,2). (3.5) 

Xi + p 

Hence, H is separable in the chart (A,/i). Moreover, H is separable also in the chart {x,y) 
and the corresponding Hamilton- J acobi equation 

H{xi, X2, dW/dxu dW/dx2) = h (3.6) 

has the complete solution W{xi,X2; h, k) = Wi{xi; h, k) + W2{x2; h, k) , Wi and W2 ful- 
filling the Sklyanin separation equations 

/l(Xi, WiiXi)) = Xih - k, /2(X2, W!,{X2)) =X2h-k , (3.7) 

with (3h + k = h . 

Proof. It is straightforward to check that the map $ : [x, y) ^ (A, fi) is a Darboux map 
for Qq] moreover, since X1—X2 = — (Ai — A2) /A1A2, the Hamiltonian H takes the canonical 



form (3.3): 



h(^x{\, /i), y{X, /i)) = (3h{x{\, /i), ?/(A, /i)) + k{x{X, fi),y{\, /x)) = (3.8) 

Al/l(^ - A -A?/il) + A2/2(^ - /?, -Xlfi2) 



Ai — A2 \ Ai A2 

/l(Ai,/ii) - /2(A2,/i2) 



Ai-A 
where 

/i(Ai, /ii) = -AiA(f - /3, -A?^i) , /2(A2, /i2) = -A2/2(f - f3, -Xlfi2) ■ (3.9) 

Ai A2 



On account of Corollary the vector field Y = Q^dH admits the qH formulation 
Y = 1/p QidK and H is separable. 

Obviously enough, H is separable also in the chart {x,y), since the map $ is a separated 
map 0], i.e., it maps separated coordinates into separated ones. Indeed, taking into 
account the form ( p.4|) of the function H, it is easily checked that the Hamilton- Jacobi 
equation H{x,dW/dx) = h has a complete solution W{xi,X2; h, k) = Wi{xi;h,k) + 
1^2(3^2; h, k), with (3 h + k = h , and that Wi, W2 fulfil the Sklyanin separation equations 
( |3.7D for the Hamilton- Jacobi equations H{x,dW/dx) = h, K{x,dW/dx) = k. o 
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4 The reduction of the tri-Hamiltonian structure of 
the Lagrange top 

If a vector field X on a manifold M is hH w.r.t. a pair of degenerate Poisson tensors 
(Po)-Pi)) a preliminary step in analysing its integrability is trying to reduce the vector 
field, its Hamiltonian functions and the Poisson tensors on a lower-dimensional manifold 
M', where one of the two Poisson tensors, say Pq? be invertible. A natural way to do 
that is to fix the values of the Casimir functions of Pq. Of course, both Pq and X can 
be properly restricted to a symplectic leaf 5*0, giving rise to a Poisson tensor Pq and to a 
vector field X' = Pq dH', H' being the restriction to 5*0 of the original Hamiltonian H. 
However, without additional assumptions. Pi is not assured to restrict to S'o, so that X' 
loses the original hH formulation. 

This situation occurs also for the tri-Hamiltonian structure of LT . Each one of the three 
Poisson tensors Pq, has two independent Casimir functions, and the generic symplectic 
leaves are four-dimensional submanifolds of M. On account of Eq.( |2.(j| ), they are 
defined as 

^0 = {m G M I = ^, ul + ul + cul - = 2a2} , (4.1) 
^1 = {m G M I = ^, ll + ll + ll = 04} , 

S2 = {m e M \ u;i7i + ^272 + c^sTs = -^«3 , li+ll + ll = 04} , 

where Oi, 02, 03 and 04 are arbitrary constants. Each Poisson tensor P^ can be properly 
restricted to a corresponding symplectic leaf S'o,, but the other two tensors do not restrict 
to the same leaf. 

Nevertheless, a quite general reduction technique given by the Marsden-Ratiu theorem 



| T9|] can be apphed; it will enables us to construct on Sa a Poisson-Nijhenuis structure 
10] induced by the tri-Hamiltonian structure on M, and on Sq a qhH formulation for 
the vector field X^. Essentially, one considers a Poisson manifold {M,P), a submanifold 
S ^ M and a distribution D C TM\g such that E := DdTS is a regular foliation with a 
good quotient J\f = S/E. Then, the theorem states that the Poisson tensor P is reducible 
to JV if the following conditions hold: 

i) the functions on M which are invariant along D form a Poisson subalgebra of C°°(M); 

ii) P{D°) CTS + D {D° being the annihilator of D in T*M). 

Analogously to previous applications of this procedure to bH structures P, ||, let us 



choose as the submanifold S a generic symplectic leaf Sa of the Poisson tensor P^ and a 
distribution D such that at each point Sa G Sa the following decomposition holds: 

T,^M = T,^Sp®D,^ , (4.2) 

Sfs being the symplectic leaf of Pp (/5 = 0, 1, 2) passing through Sq. 

This assumption assures that ii) is trivially fulfilled and that E = 0, so that the reduction 
procedure becomes a submersion H : M — > S'q onto the manifold Sa', then, it allows us 
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to endow Sa with a non degenerate tri-Hamiltonian structure, since the kernels of the 
reduced Poisson tensors P'^ vanish. Indeed, if 11* denotes the (injective) pull-back of the 
submersion 11, it is 

Ker.^P'p = iUT'{lmsn* n P^\D,^nT,^Sp)) (4.3) 

^ (n*)-i(/m,^n* n Ker.^Pp) = 0, 
where we have taken into account that 

Im^^U* C D° n Ker.^Pp = D° n (Jm.^P/,)" = D'' n (T,„5^)° ^ = ^ 0. (4.4) 

In the LT case, the distribution is as follows. 

4.1 Lemma. Let D be the distribution given by the vector fields 

d d d d 

^1 = -'^IT- + TT^' ^2 = - ^ (4.5) 
0UJ2 oujs (772 a73 

(i = a/— 1 ). Moreover, let (fi, !f2 be two generic functions. Then, for each Poisson tensor 
Pa there are two vector fields Wia and (depending on ipi and j such that 

L^,Z,+^,Z, {Pa) = ^1 A + ^2 A W2a (4.6) 

(Lz and A denoting the Lie derivative along the flow of the vector field Z and the exterior 
product of vector fields, respectively). 

Proof. It is easy to check that Lz^Pa = ZiAYija+Z2AY2ja {j = 1,2), with suitable vector 
fields Yija,Y2ja- This result, together with the identity Lfx{P) = fLx{P) + X A Pdf, 
implies ( |4.6| ), the vector fields Wja being Wja = fiYjia + ^2Yj2a + Padipj. o 

Eq. ([4.6|) implies the assumption i), since if / and g are invariant functions along D and 

Z E D, then L^z{f, g} =< df, Ly:,ziP)dg > — for each function ip. Moreover, condition 
( [4.2|) is generically satisfied as it can be easily verified. Hence, conditions i), ii) are fulfilled 
and the Marsden-Ratiu reduction technique can be applied on each symplectic leaf Sa- 
in conclusion, we have proved the following. 

4.2 Proposition. The tri-Hamiltonian structure Pp is reducible to a non degenerate 
tri-Hamiltonian structure P'^ on each one of the symplectic leaves Sa- o 

To express the reduced tensors in a particularly simple and useful form, it is convenient to 
adapt the coordinates on M to the distribution D, introducing a parametrisation including 
coordinate functions which span the subalgebra of the functions invariant along D. Let 
us choose the chart {u, v, w), related to {u, 7) by the map : M — > M : (w, 7) 1— > (u, v, w) 

ui = CUJ3- iuj2 , U2 = i'j2 - 73 , (4.7) 
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Vi = Ui, V2 = -7i, Wi = iuJ2 + CCJs, W2 = -^72 - 73 



Taking into account the tri-Hamiltonian structure Pa given by ( ^.31 ) and the definition 
( [4.1|) of 5*0, a straightforward (though lengthy) calculation allows one to verify that the 
chart {u,v) gives a parametrisation on each one of the symplectic leaves Sa', the reduced 
Poisson tensors and the tensor take the form 



P' 



( 



V - 



-1 



3 1 \ 

3 1 Ml 

-10 

Ml / 



1 








1 





\ 













-U2 






-1 













\ 





U2 








/ 



(4. 



P' 

-'2 



/ 



Ml M2 

\ M2 



-Ml 
-M2 






-U2 \ 




/ 



V ■- 



4.3 Remark. By a direct inspection, one easily concludes that the tensor N' 
(with vanishing Nijenhuis torsion) is such that P[ = N'Pq and = N'P[. 
The matrix representation of Pq and of the adjoint tensor N'* of N' are formed by Hankel 
and Frobenius blocks, respectively, so that (m, v) are Hankel- Frobenius coordinates, in the 
terminology of J^. o 



4.4 Proposition. Let us consider the map : Sa Sa '■ (m, v) (x, y) 



1 



xi = -(-Ml - yMf - 4m2) , a;2 = -(-Ml + yMf - 4m2) , (4.9) 

yi = ^(2m2 - MiMi - Vi^ul- 4m2) , 1/2 = ^(2M2 - MiMi + Mi ^M? - 4m2) • 

The chart {x,y) is a Darboux- Nijenhuis chart for the tri-Hamiltonian structure on S^; 
the reduced Poisson tensors P^ have the matrix block form 

where X = diag{xi,X2) ■ 

Proof. A straightforward computation, taking into account Eq.s ( [4 .81 ) and ( [4 .91 ). o 

(To be more precise, in order to have the Darboux-Nijenhuis chart defined in Section ^ 
one should eliminate the factor (— i) in Eg .( [4.101) , via the map x ^ ix, y ^ y). 
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5 The reduction of the vector field and the Hamilto- 
nians of the Lagrange top. 

Having established the projection of the tri-Hamiltonian structure on each one of the 
symplectic leaves Sa, the next step is to consider the reduction of the vector field and 
of the corresponding Hamiltonian functions ha- 

Unfortunately, they do not project onto Sa, since does not preserve the distribution D 
and the Hamiltonians ha are not invariant along D; hence, the tri-Hamiltonian formulation 
of Xl is lost on Sa- Nevertheless, each pair {Xl, ha) can be restrictedto the corresponding 
symplectic leaf Sa, so that Eq.f ^j]) , restricted to Sa, keeps a Hamiltonian formulation. 
Furthermore, if we consider the reduction on a symplectic leaf 5*0, we can recover, as 
a reminder of the original tri-Hamiltonian formulation, a qbH formulation for Xi; this 
suffices to provide a set of separation variables. Indeed, the following holds. 

5.1 Proposition. The vector field X^, restricted to Sq, takes the form 

Xl = PodH = -i Qo dH . (5.1) 
Its Hamiltonian H = ho^s^ takes the form 

H{x,y) = aa^H{x,y) + k{x,y) , (5.2) 

where 



H{x,y) = {f{xi,yi) - f{x2,y2)] , (5.3) 

Xl — X2 



K{x,y) = — ^ — (x2f{xi,yi) -Xif{x2,y2)) 

X^ — Xr> \ / 



X2 

m,v) = -\v' + \e+\aie+{a2+a^^)e . 

Proof. A straightforward computation. o 
On account of this result, we are just in the situation considered in Proposition with 

= crai, A = /2 = / . (5.4) 

So, Xl admits a qbH formulation; the Darboux-Nijenhuis coordinates (A, /i) are obtained 
from {x,y) via the map ( |3.5|) : 



Xi = {xi + oai) \ fJ'i = -yi {xi + craif {i = 1,2). (5.5) 



As it follows from the general results of Propositions |3.1| , |3.3| , H and K are separable 
both in the Darboux-Nijenhuis chart (A,yu) and in the chart {x,y). Using the latter, let 
us compute a solution W of the Hamilton- Jacobi equations for H and K 

dW dW dW dW 

H{xi,X2, -T—) = h, K{xi, X2, ^) = k; (5.6) 
OXi 0x2 0x1 0x2 
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taking into account the expression (|5.3|) of / and the fact that the qH function K given 



by (|373|) turns out to be = H, we have 

+ J ^/M , (5.7) 
g{0 = + are + (2a2 + ct^)^^ _ + 2(/i - aa,k). 

6 Concluding remarks 

The first result in this paper is that, reducing a la Marsden-Ratiu the tri-Hamiltonian 
structure (Pq; -Pl, P2) of onto a generic symplectic leaf S'q of each Poisson tensor, a non 
degenerate Poisson- Nijenhuis structure is obtained. The reduction depends essentially 
on the distribution D fulfilling ([4.2|) and ([4.6|) ; since D may be not unique, possibly 
different Poisson-Nijenhuis structures can be constructed on the symplectic leaf. This 
point deserves further investigations. 

The second step of the reduction procedure is the restriction of the LT vector field and 
Hamiltonian functions to the invariant submanifold 5*0, discussed in Section 5. This 
produces a quasi-bi-Hamiltonian formulation for the LT vector field and consequently, 
as a necessary outcome, a set of separation variables. An open question is whether the 
restriction of the LT vector field to other invariant submanifolds, such as the symplectic 
leaves and 5*2 of the Poisson tensors Pi and P2, gives rise to different sets of separation 
variables. 

As a last remark, we observe that the tri-Hamiltonian structure of LT has a deformation 



in the original phase space M (see, e.g., [0). In fact, there is a vector field r such that 
Lt{P2) = 2Pi, Lt-(Pi) = Pq, Lt-{Pq) = 0; in the chart (07,7) chosen in this paper, r is 
given by r = (0,0, — 2/c, cji, a;2, co^s)^. On the contrary, a recursion operator N relating 
the Poisson tensors does not exist in M. Once the reduction onto the submanifold 5*0 has 
been performed under the submersion 11 : M — 5*0, the deformation process is preserved 
since the vector field r is project able onto So- Hence, the previous relations hold for 
Pq, pi, P2 w.r.t. the projected vector field r', given by r' = —(2, ui, 0, Vi, )^ in the chart 
{u,v). As observed in Remark [4.3| , on Sq there is also an invertible recursion operator 
N' such that P{ = N'Pq and Pg = A^'P{ (consequently, one has a whole sequence of 
compatible Poisson tensors Pj+i = ^'Pj for each integer j ) . One may wonder whether 
the recursion scheme based on A^' could be inferred from the existence of the deformation 
scheme on the initial phase space, and under which conditions on the deformation vector 
field T . At the best of our knowledge, this question (which is not peculiar of LT only) has 
not yet received a satisfactory answer; in our opinion, it deserves further investigations in 
the general framework of the reduction theory for multi-Hamiltonian manifolds. 
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